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Abstract. 

I I The apparent accelerating expansion of the Universe is forcing us to examine 

r'~') the foundational aspects of the standard model of cosmology — in particular, the 

fact that dark energy is a direct consequence of the homogeneity assumption. 
We discuss the foundations of the assumption of spatial homogeneity, in the 
case when the Copernican Principle is adopted. We present results that show 
how (almost-) homogeneity follows from (almost-) isotropy of various observables. 
The analysis requires the fully nonlinear field equations - i.e., it is not possible 
Q to use second- or higher-order perturbation theory, since one cannot assume a 

^ homogeneous and isotropic background. Then we consider what happens if the 

■*^ Copernican Principle is abandoned in our Hubble volume. The simplest models 

«^ are inhomogeneous but spherically symmetric universes which do not require dark 

I I energy to fit the distance modulus. Key problems in these models are to compute 

the CMB anisotropics and the features of large-scale structure. We review how to 
^\| construct perturbation theory on a non-homogeneous cosmological background, 

^ and discuss the complexities that arise in using this to determine the growth of 

Ijr^ large-scale structure. 



o 
o 



PACS numbers: 98.80.-k, 95.36.+X, 98.80.Jk 



•£} 1. INTRODUCTION 

o 



The standard model of the Universe - tlie LCDM "concordance" model - is a 

. . perturbed FLRW model containing cold dark matter and dark energy in the form 

^^ of a cosmological constant A. This model is highly successful in being able, up to 

k> now, to fit all cosmological observations, with the same small set of parameters (^. 

%^ However, there is as yet no satisfactory explanation for the value of A, which appears 

^ in Einstein's field equations in precisely the same form as the total vacuum energy of 

quantum fields. This term is responsible for the late-time acceleration of the expansion 

of the Universe within the spatially homogeneous FLRW framework. 

The unresolved nature of A and of alternative forms of dark energy throws 
into sharp focus the foundations of the standard model - in particular, the spatial 
homogeneity assumption. We will address two questions in this paper: 

(Ql) What is the basis for the spatial homogeneity assumption? 

J Permanent address 
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(Q2) If we drop this assumption, can we find a model that is consistent with all 
cosmological observations? 

In summary, the current answer to (Ql) is that there are various results which 
motivate homogeneity (assuming the Copernican Principle); the most important one is 
that the observed ( almost- )isotropy of the CMB provides strong evidence for (almost- 
)homogeneity of the Universe on large scales [2J. But crucial open questions remain, 
which we identify and discuss. Standard cosmological perturbation theory is powerless 
here - the argument depends on a fully nonlinear analysis, since one cannot assume 
an FLRW background. 

The current answer to (Q2) is that the simplest inhomogeneous cosmological 
models - isotropic LTB (Lemaitre-Tolman-Bondi) models without any dark energy, 
and with our galaxy at the centre of a large void - are able to fit current supernova 
and some other data [3}|4l]. In this framework, the late-time acceleration becomes 
a mis-interpretation of an effect due to nonlinear inhomogeneity and curvature. In 
order to test these models against the full range of cosmological observations, it is 
necessary to develop perturbation theory on an LTB background. Although standard 
cosmological perturbation theory provides some qualitative guidance, perturbations 
on an inhomogeneous background are fundamentally different - in particular, there is 
no longer a simple separation of scalar, vector and tensor modes at first order. We 
review recent developments on this [32| , and discuss remaining open questions. 

2. WHAT IS THE BASIS FOR SPATIAL HOMOGENEITY? 

Isotropy is directly observable in principle, and indeed we have excellent data to show 
that the CMB is isotropic about us to within one part in ~ 10^ (once the dipole 
is interpreted as due to our motion relative to the cosmic frame, and removed by a 
boost). Observations of the galaxy distribution do not have the same precision, but 
there is no evidence for anisotropy. However, spatial homogeneity cannot be directly 
observed - for the simple reason that we are effectively unable to move away from 
here-and-now in order to probe spatial variations on constant time slices. The nature 
of cosmic time and length scales means that effectively our observations only access 
the past lightcone of here-and-now. Direct observation therefore cannot distinguish 
between a homogeneous distribution of matter that evolves with time down the past 
lightcone, and inhomogeneity with a different time evolution - since the past light 
cone only accesses a 2-sphere in each constant-time slice. 

The best that we can do currently is to assume a Copernican Principle - i.e., that 
our galaxy does not have an atypical view of the Universe. If we adopt the Copernican 
assumption - so that all galactic observers are assumed to see a near-isotropic CMB - 
then a nonlinear analysis gives evidence for near-homogeneity on cosmological scales, 
as we describe below. The fully nonlinear field and Boltzmann equations in covariant 
Lagrangian form are summarized for convenience in Appendix A. 

We first discuss the exact results for the case of an exactly isotropic CMB. Then 
we consider the almost-isotropic case. Finally we look at what can be deduced if 
observations of galaxies down the past light cone are isotropic. 

2.1. Exactly isotropic CMB 

What can we say if the CMB is exactly isotropic for fundamental observers? The 
pioneering result is due to Ehlers, Geren and Sachs (1968) [42]. The original EGS 
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work assumed that the only source of the gravitational field was the radiation, i.e., 
they neglected matter (and they had A = 0). We generalise their result to include 
self-gravitating matter and dark energy (extending [2|[43}|47] ) : 



(EGS) Radiation isotropy — > FLRW 

In a region, if 

• collisionless radiation is exactly isotropic, 

• the radiation four-velocity is geodesic and expanding, 

• there is dust matter and dark energy in the form of A, quintessence or 
a perfect fluid, 

then the metric is FLRW in that region. 

The fundamental 4-velocity u'^ is the radiation 4-velocity, which has zero 4- 
acceleration and positive expansion: 

.4^ = 0, e>o. (1) 

Isotropy of the radiation distribution about u^ means that (individual) photon peculiar 
velocities are isotropic for comoving observers; thus in momentum space, the photon 
distribution depends only on components of the 4-momentum p'^ along u'', i.e., on the 
photon energy E = —u^p^^: 

f{x,p) = F{x,E), F,,,...^,=OioT£>l. (2) 

In other words, all covariant multipoles of the distribution function beyond the 



monopole, defined in Eq. (A. 35 1, must vanish. In particular, as follows from 
Eq. (A.38), the momentum density (from the dipole) and anisotropic stress (from 
the quadrupole) vanish: 

qi^^O^ <" . (3) 

Equation (pi) also implies that the radiation brightness octopole Hf^^a and 
hexadecapole Ii^i,a/3 are zero. These are source terms in the anisotropic stress 



evolution equation, which is the £ — 2 case of Eq. (A. 43 1. In general fully nonlinear 
form, the tt^'^ evolution equation is 

7T<^"> + ^e^'' + ^Pra^"" + |v<^<> + 2A^^q';^ - 2w"£„^(^<)'5 
3 15 5 

+ ^aj^<)" + §v„n^''" - ll^^^n^-^ ^ 0. (4) 

Isotropy removes all terms on the left except the third, and thus enforces a shear-free 
expansion of the fundamental congruence: 

a^. = . (5) 

We can also show that u'^ is irrotational as follows. Together with Eq. (fTl), 



momentum conservation for radiation, i.e., Eq. (A. 33 1 with I — r, reduces to 

V^p, = 0. (6) 

Thus the radiation density is homogeneous relative to fundamental observers. Now we 
invoke the exact nonlinear identity for the covariant curl of the gradient, Eq. ( |A.7[ ): 

cnrW ^Pr = -2prUj^ => 6^,-0;^ = , (7) 



where we have used the energy conservation equation (A.32) for radiation. By 
assumption > 0, and hence we deduce that the vorticity must vanish: 

Wm = • (8) 
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Then we see from the curl shear constraint equation (A. 21 1 that the magnetic Weyl 
tensor must vanish: 



H, 



Eq. 



0. 



(9) 



^ actually tells us that the expansion must also be 



Furthermore, 
homogeneous. From the radiation energy conservation equation (A. 32 1, and using 
Eq. ([3]), we have Q = —'ipj./^Pr- On taking a covariant spatial gradient and using the 
commutation relation Eq. (A.8I, we find 

v^e = o. (10) 



Then the shear divergence constraint, Eq. (A. 20 1, enforces the vanishing of the total 
momentum density in the fundamental frame. 



q^^Y.1'-^ -^ E7?(/'/ 



-p/X 



(11) 



The second equality follows from Eq. ( A.29 ) , using the fact that the baryons, CDM and 



dark energy (in the form of quintessence or a perfect fluid) have vanishing momentum 
density and anisotropic stress in their own frames, i.e., 

qr = Q= ^r , (12) 

where the asterisk denotes the intrinsic quantity (see Appendix A). If we include 
other species, such as neutrinos, then the same assumption applies to them. Except 



in artificial situations, it follows from Eq. (11) that 

< = 0, (13) 

i.e., the bulk peculiar velocities of matter and dark energy [and any other self- 
gravitating species satisfying Eq. (12|] are forced to vanish - all species must be 
comoving with the radiation. 



The comoving condition ( 13 ) then imposes the vanishing of the total anisotropic 



stress in the fundamental frame: 



(l^n,") 



I 



= y^ifip}+p*jyM' = o, 



I 



(14) 



where we used Eqs. ( |A.30[ ), ( [l2| ) and p^ . Then the shear evolution equation (A. 16 1 
leads to a vanishing electric Weyl tensor 



E, 



Equations (11) 



and 



0. 



(15) 

conservation 
equation (|A.14[) and the i?-divergence constraint (|A.22[) , to homogeneous total density 



( 14 ) , now lead via the total momentum 



and pressure: 







^p = U = V,,p. (16) 

Equations (fTl), ^, ^, (10 1, (11), (14) and (16) constitute a covariant character- 
ization of an FLRW spacetime. This establishes the EGS result, generalized from the 
original to include self-gravitating matter and dark energy. (We have also provided 
an alternative, 1-1-3 covariant, analysis.) It is straightforward to include other species 
such as neutrinos. The critical assumption needed for all species is the vanishing of the 
intrinsic momentum density and anisotropic stress, i.e., Eq. (12). Equivalently, the 



energy-momentum tensor for the /-component should have perfect fluid form in the 
/-frame (we rule out a special case that allows unphysical total anisotropic stress [46]). 
The isotropy of the radiation and the geodesic nature of its 4- velocity then enforce the 
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vanishing of (bulk) peculiar velocities Vj. We emphasize that one does not need to 
assume that the matter or other species are comoving with the radiation - it follows 
from the assumptions on the radiation. 

The original EGS result (i.e., without matter or dark energy) was generalized by 
Ellis, Treciokas and Matravers (1983) fis' to a much weaker assumption on the photon 
distribution: only the dipole, quadrupole and octopole need be assumed zero. The 
key step is to show that the shear vanishes, without having zero hexadecapole - the 
quadrupole evolution equation Q no longer automatically gives afj,i, = 0, and we need 
to find another way to show this. The ETM trick is to return to the Liouvillc multipole 



equation (A.42). The £ = 2 multipole of this equation, with F^ = F^^^ = F^ua = 0, 
gives 

g|^(i^V-F_,)+i.^|^.„, = 0. (17) 

We integrate over E from to oo, and use the convergence property E^F^^^p — >• as 
i? — >■ 00. This gives 

f°° BF 
<Jc.p j E'—dE = 0. (18) 

Integrating by parts, the integral reduces to —5 /g E'^FdE. Since F > 0, the integral 
is strictly negative, and thus we arrive at vanishing shear, cr^,^ = 0. Then our proof 
above proceeds as before. 

Thus we can present a stronger generalization of the EGS result: 

(EGS-ETM) Radiation partial isotropy -^ FLRW 

In a region, if 

• collisionless radiation has vanishing dipole, quadrupole and octopole, 

^ fi — ^ flu — ^ fiva — U , V "^j 

• the radiation four-velocity is geodesic and expanding, 

• there is dust matter and dark energy in the form of A, quintessence or 
a perfect fluid, 

then the metric is FLRW in that region. 

In summary, the EGS-ETM result, suitably generalized to include baryons and 
CDM and dark energy, is the most powerful basis that we have - within the framework 
of the Copernican Principle - for background spatial homogeneity and thus an FLRW 
background model. 

Although this result applies only to the "background Universe", its proof 
nevertheless requires a fully nonperturbative analysis. 

2.2. The real universe: almost-isotropic CMB 

In practice we can only observe approximate isotropy. Is the EGS result stable - i.e., 
does almost-isotropy of the CMB lead to an almost-FLRW Universe? This would 
be the realistic basis for a spatially homogeneous Universe (assuming the Copernican 
Principle). It was shown to be the case, subject to further assumptions, by Stoeger, 
Maartens and Ellis [2l|49l: 
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(Almost-EGS) Realistic basis for spatial homogeneity 

In a region of an expanding Universe with cosmological constant, if all 
observers comoving with the matter measure an almost isotropic distribution 
of collisionless radiation, and if some of the time and spatial derivatives of the 
covariant multipoles are also small, then the region is almost homogeneous. 

We emphasize that the perturbative assumptions are purely on the photon 
distribution, not the matter or the metric - and one has to prove that the matter 
and metric are then perturbatively close to FLRW. Once again, a nonperturbative 
analysis is essential, since we are trying to prove an almost-FLRW spacetime, and we 
cannot assume it a priori. 

Almost-isotropy of the photon distribution means that 

F^^...^^ix,E)=0{e), £>1, (20) 

where e is a (dimensionless) smallness parameter. The brightness multipoles Umi 
have dimensions of energy density and we therefore normalize them to the monopole 
n = pr/4:Tr: 

^=Oie). (21) 

The task is to show that the relevant kinematical, dynamical and curvature 
quantities, suitably non-dimensionalized, are 0{e). For example, the dimensionful 
kinematical quantities may be normalized by the expansion, cr^^/OjWp/G. The proof 
then follows the same pattern as our proof above of the exact EGS result - except 
that at each stage, we need to show that quantities are C(e) rather than equal zero. 

However, in order to show this, we need smallness not just of the multipoles, 
but also of some of their derivatives. Smallness of the multipoles does not directly 



imply smallness of their derivatives, and we have to assume this 47 50 . It has been 



claimed 47 , 50 that the derivatives are generically not small, but this has not been 
rigorously established - which is not surprising, since it is a difficult problem and may 
be related to the averaging problem in cosmology. The issue needs to be settled by 
further analysis. 

If all observers measure small multipoles, then it may be possible to show 
- perhaps using other observations - that the time and space derivatives on 
cosmologically significant scales must also be small. A number of experiments has 
been proposed to test the Copernican Principle by looking for violations of isotropy at 
events down our past lightcone. These include looking for spectral distortions of CMB 
photons scattered by ionized gas [5H[52] : such distortions are induced by anisotropics 
in the CMB as seen by distant observers, and so provide in principle a neat way of 
confirming the Copernican assumption as used here. A similar test uses the kinematic 



Sunyaev-Zeldovich effect in clusters to observe the dipole around distant observers 27 
The almost-EGS result then gives a framework for probing inhomogeneities via such 
observations. Indeed, these tests may provide a way of constraining spatial gradients 



of the low-^ multipoles 54 



In addition, it may be possible to strengthen the almost-EGS result above by 
proving that it is sufficient for only the first 3 multipoles and their derivatives to be 
small. This would be an almost-EGS-ETM result, and would represent a more realistic 
foundation for almost-homogeneity than the almost-EGS result. 

Note also that the almost-EGS result has not been proven for quintessence or 
perfect fluid dark energy, and this needs further investigation [54 . 
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2.3. Isotropic matter distribution on the past lightcone 

We have seen the power of an (almost-) isotropic CMB as the fundamental basis for 
a spatially homogeneous background model of the Universe. What can we say if the 
matter distribution is (almost-) isotropic on the past lightcone of the observer? 

A "matter" counterpart to the EGS result is based on analysis of the observations 
of galaxies down the past lightcone of the observer. Without assuming the Copernican 



Principle, we have the following isotropy result 55 - 57 

Matter lightcone- isotropy — > spatial isotropy 

If one fundamental observer comoving with the matter measures isotropic 
area distances, number counts, bulk peculiar velocities, and lensing, in an 
expanding dust Universe with A, then the spacetime is isotropic about the 
observer's worldline. 

The original result neglected A and CDM, and we can incorporate both in the 
same way as in the EGS case. Note that isotropy of (bulk) peculiar velocities seen 
by the observer is equivalent to vanishing proper motions (tranverse velocities) on the 
observer's sky. Isotropy of lensing means that there is no distortion of images, only 
magnification. 

The proof of this result requires a non-perturbative approach - there is no 
background to perturb around. Since the data is given on the past lightcone of the 
observer, the covariant Lagrangian approach used in the previous subsection is not 
suitable. Instead, we need the fully general metric, adapted to the past lightcones 
of the observer worldHne C. We define observational coordinates x^ = [w,y,9,(p), 
where x"" = {0, <j>) are the celestial coordinates, w = const are the past light cones on C 
(y = 0), normalized so that w measures proper time along C, and y measures distance 
down the light rays {w,6,(j)) = const. A convenient choice for y is y = 2; (redshift) 
on the lightcone of here-and-now, w = wq, and then keep y comoving with matter 
off the initial lightcone, so that u^ = 0. Then the matter 4-velocity and the photon 
wave-vector are 

u^ = (l + z) (1,0,1^'^), k^ = w^^, l + z = u^k^, (22) 

where V"" — dx"" /dw are the transverse velocity components on the observer's sky. 
The metric is 

ds^ = - A^dw^ + 2Bdwdy + 2Cadx''dw + D^{dn^ + L„fcda;Ma;'')(23) 

A^ ={l + z)-^ + 2CaV^+g,bV''V\ B^^, (24) 

dy 

where D is the area distance, and Lab determines the lensing distortion of images via 
the shear of lightrays, 

^'''=2B^f (2^) 

The number of galaxies in a solid angle dO and a null distance increment dy is 

dN = Sn{l + z)D^Bdndy, (26) 

where S is the selection function and n is the number density. 

Before specializing to isotropic observations, we identify how the observations in 
general and in principle determine the geometry of the past light cone w — Wq oi 
here-and-now, where y — z: 
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• Area distances directly determine the metric function Z?(wo,z,x°). 

• The number counts (given a knowledge of S) determine Bn and thus, assuming a 
knowledge of the bias, they determine /i,„(K;o, z, a;°) = B{wo, z, x"')pmiwo, z, x""), 
where pm = Pb + Pc is the total matter density. 

• Transverse (proper) motions in principle directly determine V°'{wo, z,x^). 

• Image distortion determines Lab(wo, ZjX"^). (The differential lensing matrix aab 
is determined by Lab, D, B.) 



Then 55 56 



Lightcone observations —>■ spacetime metric 

Observations {D,N,V"',Lab) on the past lightcone w — wq determine 
in principle {gab,u'^, Pm = Bpm) on the lightcone. This is exactly the 
information needed for Einstein's equations to determine _B, CaOnw = Wq, so 
that the metric and matter are fully determined on the lightcone. Finally, the 
past Cauchy development of this data determines .gp,y,u^,p„i in the interior 
of the past lightcone. 

If we assume that observations are isotropic, then 

Momentum conservation and the yy field equation then give the following equations 
on w = Wq (55JI57]: 

C, = (1 + zy' / (1 + z)B^adz (28) 

Jo 



dz 



(29) 



(1 + z) Dpradz 
'0 

where a prime denotes d/dz. These imply that Ba ~ ^ Ca, so that pjn,a = - 
and hence the metric and matter are isotropic on w ^ wq. This can only happen if 
the interior of w — wq is isotropic. If observations remain isotropic along C, then the 
spacetime is isotropic. 

Adopting the Copernican Principle, it follows that all observers see isotropy and 
therefore that spacetime is isotropic about all galactic worldlines - and hence spacetime 
is FLRW. This result then becomes a "matter" alternative to the EGS, as a basis for 
FLRW: 

Matter lightcone- is otopy -^ FLRW 

In an expanding dust region with A, if all fundamental observers measure 
isotropic area distances, number counts, bulk peculiar velocities, and lensing, 
then the spacetime is FLRW. 

In essence, this is the Cosmological Principle, but derived from observed isotropy 
and not from assumed spatial isotropy. Note the significant number of observable 
quantities required. 

We can actually give a much stronger statement than this, based only on distance 
data. An important though under-recognized theorem due to Hasse and Perlick tells 



us that 58 
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(HP) Isotropic distances -^ FLRW 

If all fundamental observers in an expanding spacetime region measure 
isotropic area distances up to third-order in a redshift series expansion, then 
the spacetime is FLRW in that region. 

The proof of this relies on performing a series expansion of the distance-redshift 
relation in a general spacetime, using the method of Kristian and Sachs [59 , and 
looking at the spherical harmonic multipoles order by order. We illustrate the proof 
in the 1+3 covariant approach (in the case of zero vorticity). Performing a series 



expansion in redshift of the distance modulus, we have, in the notation of 60 
m- M 



25 



5 logio ; 
5 



log 






10' 



{K^KPV^Up 



J o 



R^.K^K'' 



^{KPK'nVpUnY 



A+0{z'), (30) 



where 



^M 



kP 



'k 



and A'^l 



-{u^' 



\o ' 



(31) 



v\o 



denotes a past-pointing null vector at the observer O in the direction e^. Comparing 
with the standard FLRW series expansion evaluated today, we define an observational 
Hubble rate and deceleration parameter as 



H' 



obs I 



= Rf^K^Vn 



obs I 

y lo 



(32) 
(33) 



(K-KPV^Up)^ 

We can also give an effective observed cosmological constant parameter from the 0{z'^) 
term: 

5 ,. „K.x . RpuK^'K'' 



nt' = 



[l-q^^^)-5 



12(i?°bs)2 



6(i^obs)3 



(34) 



The argument of 58 relies on proving that if all observers measure these three 



quantities to be isotropic then the spacetime is necessarily FLRW. In a general 
spacetime 



60 



HI 



obs 



^^e 



Apc^ 



' ^lU'^ 



(35) 



where A^e'' is a dipole and cr^jye^e'^ is a quadrupole. Hence, if all observers measure 
Hq^^ to be isotropic, then a^^ = Q = A^. In a spacetime with isotropic H^^^ the 



generalized deceleration parameter is given by f81' 



(^0 



obs\ 



90 



obs 



1 



1 



1 



-A 



1 



2?^ 



3v^e 



E. 



(m^; 



- 2^{i^-^) 



(36) 
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If the dipolc of this term vanishes then we see from Eq. (A. 20 1 that the energy flux 



must vanish as weU as spatial gradients of the expansion. Excluding models with 



unphysical anisotropic pressure, Eq. (A.17) then shows that the electric Weyl tensor 



must vanish, and it follows that the spacetime must be FLRW. The more general 



proof in 58 uses fl'^^ to show that the vorticity must necessarily vanish along with 



the anisotropic pressure. 

It is an open but important question how this result translates to the case of 
almost-isotropy. 

3. COSMOLOGY WITHOUT THE COPERNICAN PRINCIPLE 

While the Copernican Principle remains untested, inhomogeneous models should be 
explored. As we discussed in the previous section isotropic observations imply spherical 
symmetry in the presence of dust matter and A, leading to the LTB model. A more 
general spherically symmetric perfect fluid solution is the Lemaitre model [61| , which 
can include inhomogeneous radiation, or other types of barotropic perfect fluid. LTB 
solutions are the dust subcase. An interesting explanation for the dark energy problem 
in cosmology is one where the underlying geometry of the universe is significantly 
inhomogeneous on Hubble scales. Spacetimes used in this context are usually LTB 
models - so-called void models, first introduced in Is]. These models can look like dark 
energy because we have direct access only to data on our lightcone and so we cannot 
disentangle temporal evolution in the scale factor from radial variations. Void models 
are considered ungainly compared with standard cosmology because naively they 
revoke the Copernican Principle, placing us at or very near the centre of the universe. 
It may be that the models used are very simplistic descriptions of inhomogeneity, and 
more elaborate inhomogeneous ones will conform to some version of the Copernican 
Principle and still satisfy observational constraints on isotropy. 

Instead, we may think of LTB void models as smoothing all observables over the 
sky, thereby compressing all inhomogeneities into one or two radial degrees of freedom 
centred about us - and then we avoid placing ourselves 'at the centre of the universe' 
in the standard way. In this sense, these models are a natural first step in developing 
a fully inhomogeneous description of the universe. However, given the results of the 
last section it is not clear if this interpretation will work. Indeed, if we can extend the 
exact result that matter lightcone-isotopy — ^ FLRW, to the almost-isotropic case, this 
could rule out radially inhomogeneous models with this interpretation. Furthermore, 
it has not been shown that a sky-averaged inhomogeneous model would give an LTB 
dust solution. It is likely that this is not the case because averaging and evolution do 
not commute. 

Alternatively, within the multiverse context, one can imagine a vast universe in 
which our observable patch happens to be rather inhomogeneous. We can even imagine 
a multiverse in which there are many void-like regions; even if we happened to be near 
the centre of one with a Hubble-scale inhomogeneity, this may be natural within a 
larger context - perhaps in the same way we discovered that the Milky Way is not 



particularly special once understood in the context of a plethora of galaxies 62 . With 
this idea, we need not violate the Copernican Principle if we live near the centre of 
a Hubble scale void; rather, we should just change our perspective, and think of the 
Hubble scale as very small rather than very large. 



Indeed, as argued in 62 , the anthropic 'explanation' for the current value of A, 



which relies on a multiverse of some sort for its philosophical underpinning, necessitates 
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the violation of the Copernican Principle - simply because the vast majority of universe 
patches are nothing like ours, and not at all suitable for complex life. 

3.1. Hubble-scale voids 

An inhomogeneous void may be modelled as a spherically symmetric LTB model with 
metric 

ds^:^~dt^+ ^ " ] ^ dr^ + a\{t,r)r'^dn'^ , (37) 

1 — K{r)r-^ 

where the radial (ay) and angular (aj^) scale factors are related by ay = d{a±r)/dr. 
The curvature k = K{r) is not constant but is instead a free function. The two scale 
factors define two Hubble rates: 

Hj_ = Hj_{t,r) = —, Hn=Hn{t,r) = ^, (38) 

aj_ ail 

and the Friedniann equation takes on its familiar form: 

--^ = n„,al^ + n,al^, (39) 

where rJm(r) + ^^^(r) = 1 and r2m(r) is a free function, specifying the matter density 
parameter today. In general, H±g(r) is also free, but removing the decaying mode by 
enforcing a uniform bang time [63[|64| fixes this in terms of flmir). 

With one free function we can design models that give any distance modulus we 
like. If we choose il.m{i") to reproduce exactly a LCDM D{z), then the LTB model is 
a void with steep radial density profile which is, strictly speaking, non-differentiable 



at the origin if we want ^o < 13 (see also 34 ). Much has been made of this non- 
differentiability, but it is irrelevant for this cosmological modelling. In essence, the fact 
that a sharp-profiled void can reproduce the LCDM D{z) exactly is of mathematical 
interest only; a void model should be constrained directly by data, and not matched 
to a best fit LCDM model. If this is done, a smooth void is perfectly acceptable 



with present data 35 , but this provides a feasible route to distinguish between a 



smooth void and LCDM 30 . However, it is impossible to tell the difference between 



an evolving dark energy FLRW model and a void model, using distance data alone. 

3.2. Perturbation theory in void models 

An important open problem in inhomogeneous models is the modelling of structure 
formation. This is important partly because it provides a means for distinguishing 
between FLRW and LTB. One example of where we might see an effect is in the 
peak in the two-point matter correlation function attributed to the Baryon Accoustic 
Oscillations (BAO). It has been shown that if LTB perturbations evolve as in FLRW, 



then BAO can be decisive in ruling out voids [26 29 . Whether this assumption is 
valid however requires a full analysis of perturbations. 
There have been two approaches so far: 

(i) Using a covariant H-l-t-2 formalism which was developed for gauge-invariant 



perturbations of spherically symmetric spacetimes 65 66 . The full master 



equations for LTB have not yet been derived, but some progress has been made in 



the 'silent' approximation, neglecting the magnetic part of the Weyl tensor 64|67 
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(ii) Using a 2+2 covariant formalism 68 69 , developed for stellar and black hole 



physics. The full master equations for LTB perturbations were presented in 32 
(see also [TO]). 

In FLRW cosmology, perturbations split into scalar, vector and tensor modes that 
decouple from each other, and so evolve independently (to first order). Such a split 
cannot usefully be performed in the same way in a spherically symmetric spacetime, 
as the background is no longer spatially homogeneous, and modes written in this 
way couple together. Instead, there exists a decoupling of the perturbations into two 
independent sectors, called 'polar' (or even) and 'axial' (or odd), which are analogous, 
but not equivalent, to scalar and vector modes in FLRW. These are based on how the 
perturbations transform on the sphere. Roughly speaking, polar modes are 'curl' free 
on S"^ while axial modes are divergence free. Further decomposition may be made 
into spherical harmonics, so all variables are for a given spherical harmonic index 
i, and modes decouple for each £ - analogously to fc-modes evolving independently 



on an FLRW background. A full set of gauge-invariant variables were given by 69 
who showed that there exists a natural gauge - the Regge- Wheeler gauge - in which 
all perturbation variables are gauge-invariant (rather like the longitudinal gauge in 
FLRW perturbation theory). Unfortunately, the interpretation of the gauge- invariant 
variables is not straightforward in a cosmological setting. 

Most of the interesting physics happens in the polar sector, so we will discuss 
that case, following |32'. The general form of polar perturbations of the metric can 
be written, in Regge- Wheeler gauge, as 

ds^ = - [1 + (2?7 - X - (p)Y] dr " dtdr 

V 1 — Kr^ 

+ [1 + (X + ^)r] -J +alr^{l + ^Y)dn^, (40) 

(1 — Kr^) 

where r](t,r), x(^j''): ¥'(^,?') and <;{t,r) are gauge-invariant variables. The notation 
here is such that a variable times the spherical harmonic Y has a sum over i, m, 
e.g., ipY — J2'iLoJ2m=-e^^"^i^^)^i"i(^'^)^ where x° are coordinates on S*^, and 
X = (t,r). The general form of polar matter perturbations in this gauge is given 



by 



UA + 1 wflA + -hABU^ ] Y,vYa 



(41) 

p =/^^(l + Ay), (42) 

where u, w and A are gauge-invariant velocity and density perturbations and Hab 
is the metric perturbation in the x part of the metric; a colon denotes covariant 
differentiation on the 2-sphere. The unit vectors in the time and radial directions are 



\A _ n n-\ ^A 



U 



Vi^ 



2 



nr^ 



= (1,0), n''=\Q,- . (43) 



The elegance of the Regge- Wheeler gauge is that the gauge-invariant metric 
perturbations are master variables for the problem, and obey a coupled system of 
PDEs which are decoupled from the matter perturbations. The matter perturbation 
variables are then determined by the solution to this system. We outline what this 
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system looks like for ^ > 2; in this case 77 = 0. The generahzed equation for the 



gravitational potential is 32 



ip + 4:Hj_ip - 2^-ip = S^{x, <^)- 



(44) 



The left hand side of this equation has exactly the form of the usual equation for a 
curved FLRW model, except that here the curvature, scale factor and Hubble rate 
depend on r. On the right, S^ is a source term which couples this potential to 



gravitational waves, Xj and generalized vector modes, 
are sourced by ip: 



<?■ 



These latter modes in turn 



X + x''-iHnx-2Wx' + 



16nGp - 



6M 



cj + 2i/j,? = -x'. 
The prime is a radial derivative defined by 
, ^/l — Kr^ 



{i-m+2) 



a\r'^ 



X = 5'x(<r,(y5), 



(45) 
(46) 

(47) 



The gravitational field is inherently dynamic even at the linear level, which is not 
the case for a dust FLRW model with only scalar perturbations. Structure may grow 
more slowly due to the dissipation of potential energy into gravitational radiation. 
Since H± = H±{t,r), a± = a±{t,r) and k — K{r), perturbations in each shell about 
the centre grow at different rates, and it is because of this that the perturbations 
generate gravitational waves and vector modes. This leads to a very complicated set 
of coupled PDEs to solve for each harmonic i. 

In fact, things are even more complicated than they first seem. Since the scalar- 
vector-tensor decomposition does not exist in non-FLRW models, the interpretation 
of the gauge-invariant LTB perturbation variables is subtle. For example, when we 
take the FLRW limit we find that 



ip= -2'^-2HV-2 



(1 - Kr2) 



hr 



1 



/l(^) 



-ndr 



(l-Kr^) £{£ + 1) - 4{l - Kr^) 



-a, 



2r2 



h<^^^\ (48) 



where ^ is the usual perturbation space potential, V is the radial part of the 
vector perturbation, and the /I's are invariant parts of the tensor part of the metric 
perturbation. Thus if contains scalars, vectors and tensors. A similar expression 
for <; shows that it contains both vector and tensor degrees of freedom, while x is 
a genuine gravitational wave mode, as may be seen from the characteristics of the 
equation it obeys. This mode mixing may be further seen in the gauge-invariant 
density perturbation which appears naturally in the formalism: 

SnGpA = -if" - 2Wip' + (i?|| + 2Hj_)ip + Wx' + Hj_x 

'iie + 1) 



{i~m +2) 

2a\r 



2Hi+4:H\\H^-8nGp (x + 'p) 
-X + 2HW + 2{m+H^)W<;, 



(49) 
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where 



W 



VT^ 



(50) 



a±r 

When evaluated in the FLRW Umit the mode mixing becomes more obvious stiU: A 
contains both vector and tensor modes, while its scalar part is 

AnGa^p A = V^* - 3^9^* - 3{n^ - k)*, (51) 

which giv es the usual gauge invariant density fluctuation 5p^'~^^' ^ 6p + drp{B — 



OtE) 71 . Here, V refers to the Laplacian acting on a 3-scalar, so that 



(1 - Kr^)d^ + 



(2-3«;r2)^ £{£+!) 



(52) 



(which also helps one identify a Fourier k-mode in this context). The fact that A is 
more complicated is because the gauge-invariant density perturbation includes metric 
degrees of freedom in its definition; gauge-invariant variables which are natural for 
spherical symmetry may not be natural for homogeneous backgrounds. A gauge- 
dependent A may be defined which reduces to Sp'^'^^^ in the FLRW subcase, but its 
gauge-dependence will cause problems in the inhomogeneous case. 

A gauge-invariant variable which reduces to a pure scalar in the FLRW limit 
is 



32 



C = A" + 21^A' - 



£{t + l) 



A + rWi' + r 2,W^ - 



1 



L 



where 



A 



87rGpa_L [i?I^A 



M 



(53) 



(54) 



In the FLRW limit, A contains no tensors. Furthermore, 



3a j^ 
2W 



h^r-iU^'^^-'^ 



^(£+l)-3 



a^j-"^ 



iW 



is a pure vector in the FLRW limit. In the FLRW limit, we find 

_ P 
where the curvature perturbation is 7^ = *-'H('H$ + 9r*)/(9r'H-'H^- 



c 



n 



37^ - 3* 



(55) 



(56) 



-k) [71]. These 

generalized scalar-vector-tensor variables may be useful in interpreting perturbations 
of void models. 

These equations have not yet been solved in full generality. We expect different 
structure growth in LTB models, but it is not clear what form the differences will take. 
It seems reasonable to expect that the coupling between scalars, vectors and tensors 
will lead to dissipation in the growth of large-scale structure where the curvature 
gradient is largest, as it is the curvature and density gradients that lead to mode 
coupling. In trying to use structure formation to compare FLRW to LTB models, 
some care must be taken over the primordial power spectrum and whatever early 
universe model is used to generate perturbations - since there is a degeneracy with 
the primordial power spectrum and the features in the matter power spectrum. 
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3.3. The small-scale CMB in void models 

The physics of decouphng and hnc-of-sight effects contribute differently to the 
CMB, and have different dependency on the cosmological model. In sophisticated 
inhomogeneous models both pre- and post-decoupling effects will play a role, but 
Hubble-scale void models allow an important simplification for calculating the 
moderate to high £ part of the CMB. 

The comoving scale of the voids which closely mimic the LCDM distance modulus 
are typically O(Gpc). The physical size of the sound horizon, which sets the largest 
scale seen in the pre-decoupling part of the power spectrum, is around 150 Mpc. This 
implies that in any causally connected patch of the Universe prior to decoupling, the 
density gradient is very small. Furthermore, the comoving radius of decoupling is 
larger than 10 Gpc, on which scale the gradient of the void profile is small anyway 
(by assumption). For example, at decoupling the total fractional difference in energy 



density between the centre of the void and the asymptotic region is around 10% 72 
hence, across a causal patch we expect a maximum 1% change in the energy density 
in the radial direction, and much less at the radius of the CMB that we observe for a 
Gaussian profile. This suggests that before decoupling on small scales we can model 
the universe in disconnected FLRW shells at different radii, with the shell of interest 
located at the distance where we see the CMB. This may be calculated using standard 



FLRW codes, but with the line-of-sight parts corrected for 29 33 . 

For line-of-sight effects, we need to use the full void model. These come in two 
forms. The simplest effect is via the background dynamics, which affects the area 
distance to the CMB, somewhat similar to a simple dark energy model. This is the 
important effect for the small-scale CMB. The more complicated effect is on the largest 
scales through the Integrated Sachs- Wolfe effect (see [73] for the general formulas in 
LTB). This requires the solution of the perturbation equations presented above. 



The CMB can be fit in void models in different ways. In 29 33 , it was shown 
that the CMB can be very restrictive on void models, although [33| showed that 
with a varying bang time, the data for Hq, SNIa and CMB can be simultaneously 
accommodated. Including inhomogeneous radiation in the background, it was argued 
in 72 that the CMB can be accommodated along with other local observations with 
a homogeneous bang time, but with asymptotic curvature at the CMB radius. It is an 
open question exactly what constraints the small-scale CMB places on a generic void 



solution. An example from 72 of how closely a void model can reproduce the CMB 
power spectrum found in a concordance model is shown in Fig. [Tl It is remarkable 
that the void models can reproduce the LCDM CMB power spectrum so closely. As 



emphasised in 33 , it is often taken for granted that the CMB tells us that the universe 
is close to flat; these examples show that curvature can in fact be very large, but just 
not homogeneous. 

3.4. Summary of inhomogeneous models 

An inhomogeneous LTB void model, even if it over-simplifies nonlinear inhomogeneity, 
does produce some rather remarkable results. The apparent acceleration of the 



universe can be accounted for without dark energy. It has also recently been shown 72 
that LTB models can simultaneously solve the so-called lithium problem, the apparent 
mismatch between ^Li abundance measured locally (thereby probing our past worldline 
at nucleosynthesis), and that inferred from measurements of deuterium and the 
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Figure 1. The normalized CMB angular power spectrum, from [72] . The power 
spectrum is shown against a default flat concordance model with zero tilt. There 
is nothing between the two models for high £, with the maximum difference around 

1%. 



baryon-photon ratio at the CMB (which probe nucleosynthesis at co-moving scales 
O(Gpc) away from us). In the concordance model, the lithium problem is a 2.4 — 5.3(7 



discrepancy 80 , yet it can easily be accommodated with Hubble scale inhomogeneity 



in the baryon-photon ratio at early times. 

It has also been shown that the small-scale CMB need not over-constrain voids 
and so rule them out 33 72 , but it seems that the simplest possibility - homogeneous 
bang time and asymptotic flatness - may not be possible. 

The major open problems in inhomogeneous solutions to the dark energy problem 
are: 

• CMB. The details of an inhomogeneous radiation era still have to be investigated. 
To calculate the large scale CMB, there are several additional effects which must 
be taken into account. The most important is the Sach-Wolfe effect, which 
requires knowing the perturbation spectrum at the time of decoupling on the 
largest angular scales. While it might not deviate from FLRW for a central 
observer because of the reasons discussed above, this needs further investigation. 
In addition, the Integrated Sachs- Wolfe effect will contribute to the line-of-sight 



part of the CMB calculation 73 74 . 

Structure formation. The complexity of the perturbation equations is a major 
stumbling block. Attempts have been made to solve a limited subcase of the 
equations 32]|64 67 , but it is not clear how close those approximations are to the 
full solution. This is required for a reliable estimate of corrections to the power 
spectrum and to the BAO peak. 

Inflation. How can a void model arise in the first place? What kind of voids are 
possible from an inflationary era? Is there a natural way to produce many voids, 
and so ease the Copernican problem? Standard inflation can't easily give a single 
void, but one can imagine an inflationary potential -- or perhaps a separate earlier 
period of inflation some 50 efolds before the standard slowly rolling inflationary 
phase - which imprints an extra scale in the primordial power spectrum. 

The Copernican problem. As we have discussed above, it is not easy to flnd 
an inhomogeneous model with nearly isotropic observations which satisfies the 
Copernican Principle. Further work is needed to see if it is possible to avoid or 
understand the very unsatisfactory fine-tuning problem which exists in LTB voids: 
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that we are within tens of Mpc of the centre 2T|[22 36 . One way one can imagine 



this question being addressed is some kind of anthropic selection effect which 
makes inteUigent Hfe more hkely in regions where the dark matter concentration 
is lowest JU For example, one could imagine that dark matter inhibits the formation 
of stable solar systems of sufficient longevity. Therefore, a naive calculation of a 
1 in ~ 10^ chance of us being within tens of Mpc of the centre of a Hubble-scale 
void would have to be weighted appropriately, thereby reducing this fine-tuning 
problem. 

4. DISCUSSION 

Without a theoretical understanding of the value of the cosmological constant the 
concordance model remains phenomenological no matter how strongly observations 
appear to support the model. The concordance model suffers an important coincidence 
problem: that dark energy starts to dominate at around the time our solar system 
forms. Until we demonstrate observationally that the Universe is homogeneous 
on large scales, we should consider inhomogeneous spacetimes even if they are 
philosophically uncomfortable, particularly in light of the fact that in their simplest 
incarnation they can explain away the dark energy problem through inhomogeneity, 
without apparently causing other problems. It is worth reflecting on how remarkable 
this is. Measurements of a nonzero A may in fact be thought of as an explicit 
consequence of the Copernican assumption. 

It is important therefore to re-examine the basic assumptions of the standard 
concordance model - both to test the Copernican Principle at the heart of modern 
cosmology, and to understand what freedom we have to develop new cosmological 
models that may be able to explain cosmic acceleration using curvature and 
inhomogeneity rather than a cosmological constant. In order to do this, we need 
non-perturbative methods and also new types of nonlinear perturbations where 
the nonlinearity arises from the inhomogeneous background. Moreover, the whole 
inflationary paradigm needs to be re-examined to see if inhomogeneities this large 
could arise naturally. 

We have shown how the Copernican assumption when combined with the high 
isotropy of the CMB implies FLRW under fairly weak assumptions - mainly, that 
the CMB rest frame is geodesic. Provided that the dark matter is CDM that is 
comoving with the CMB, this geodesic property will hold. But in interacting dark 
energy-dark matter models, for example, the geodesic radiation frame is an additional 
assumption. In order to test the Copernican Principle at the foundation of the FLRW 
models, observations which look inside our past lightcone to estimate how distant 
observers see the CMB 27 51 ■ 53 are necessary. In order to observationally prove 
FLRW, however, we also must demonstrate that galaxies follow geodesies, and that 
dark energy is comoving with the CMB. 

If galaxies were non-geodesic, then this would leave a signature as a dipole in 
the Hubble law which grows linearly with distance (and so is distinguishable from a 
normal bulk flow, which is constant) [45| which is derived from the term which is linear 
in redshift in the area distance-redshift relation - see Eq. ( pS] ). Possible dark energy 
flux would leave a similar signal in the quadratic term in the area distance-redshift 
relation. The generalized deceleration parameter, in a fully inhomogeneous spacetime 

§ S. Rasanen, private communication. 
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is, following 81 , 



H^o 



-p A 

2^ 3 

2 



^V,A- 



12 



-CTui^cr' 






A'a,, 



rq^L 



3 

5^- 



P (v,w^ + 2A,uj^ 



+e''e- 



-2V(^Av) 



^(p^^) 



£;, 



ifii^) 



i-(p,.) 



-'^(p^.) - 3^^(p^, 



-5A 



(/.cr^-y) 



9 



(m'^i't) 



Be^e'^e^eV 



{f.LV^'yd) • 



(57) 

If it were confirmed observationally that the CMB frame is geodesic, and that the CMB 
is almost isotropic around distant observers, then this H^q^ should be an isotropic 
observable - a potentially useful consistency check. 

Although it seems obvious that CDM and baryons are just geodesic dust, we 
have to be careful because of issues highlighted from the averaging problem. While 
a distribution of non-relativistic CDM particles obey a dust equation of state at all 
times, more and more of the CDM particles are locked up as galaxies and clusters form, 
and yet those virialized clusters then become the 'dust particles', within the standard 
paradigm. Are these geodesic? It is not yet clear how to treat this situation properly, 
and it is conceivable that there is a backreaction effect which changes the effective 
equation of state of matter at late times [82] (amongst other effects) . Depending on 
precisely what is calculated, this backreaction can be divergent in the concordance 



model 83 , and remains a significant problem to calculate non-linearly (although 
see 



84 who find no such divergences using a different method of averaging). 



Off-lightcone observations will be able to place limits on the spatial gradients of 
the CMB multipoles as well as the multipoles themselves, which we have shown is 
also a requirement for deducing almost-FLRW geometry from almost-isotropic CMB 
(the almost-EGS result). The important exact ETM result states that if the first 
three multipoles are zero then all the rest must be too. It seems reasonable that this 
would extend to an almost-ETM result if the first 3 multipoles are small. In this case 
observations which look inside our past lightcone might only need to confirm that the 
first three multipoles around distant observers are small in order to observationally 
confirm FLRW. These issues deserve further investigation, and are important because 
it might be possible to construct a genuinely inhomogeneous model which solves the 
dark energy problem, but which also satisfies the Copernican principle, maintaining 
the isotropy of the CMB for all or 'most' observers (as claimed in [851, for example). 

Dropping the Copernican Principle, we have the key exact result that isotropy 
of matter observations down our past lightcone implies an LTB geometry. However, 
it is an open question whether almost-isotropy of matter observations on the past 
lightcone implies an almost-spherically symmetric LTB model. We have discussed 
a complete analysis of structure growth at first order, which is the next step in 
developing LTB models, and shown how this is much more involved than the linear 
perturbations of FLRW. However, we argue that it is worth the effort. All tests 
of the Copernican Principle, including those which can be done by examining the 
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background dynamics down our past lightcone 86p7 , are difficult and sensitive. Only 



by developing a family of inhomogeneous spacetimes to the same level of sophistication 
as the standard concordance model, and directly comparing them side by side, will 
we really be able to understand whether A is real, or actually a consequence of our 
homogeneity assumption. 

Acknowledgments: We thank Timothy Clifton, Ruth Durrer, George Ellis, Pedro 
Ferreira, David Matravers, John Moffat, Syksy Rasanen and Marco Regis for 
comments and/or discussions. CC is supported by the National Research Foundation 
(NRF, South Africa) and RM is supported by the UK Science & Technology Facilities 
Council. This work was partially supported by the South African National Institute 
for Theoretical Physics and by a Royal Society (UK)/ NRF exchange grant between 
the Universities of Portsmouth and Cape Town. 



Inhomogeneity and the foundations of concordance cosmology 20 

Appendix A. Fully nonlinear field and Boltzmann equations 



We summarize, using 88 , the key equations in tlie 1+3 covariant Ehlers-Ellis 



formalism (see also [89] in this volume, and 90,91] for reviews). This formalism 
provides a physically transparent formulation of the field equations and the Boltzmann 
equation in full nonlinear generality. The formalism leads to a simple 'top-down' 
approach to linear perturbations around FLRW - i.e., one starts from the nonlinear 
spacetime and then takes the FLRW limit. By contrast, the 'bottom-up' approach of 



the usual coordinate-based perturbation theory (see 71 for a review), starts from 
FLRW and perturbs away from it. In the context of perturbed FLRW models, 
the covariant approach has no advantage over the coordinate-based approach. But 
beyond the FLRW framework, when one cannot assume there is an FLRW background 
to perturb away from, the coordinate-based approach cannot be applied, while the 
covariant approach is applicable. 

The Ehlers-Ellis formalism is a covariant Lagrangian approach to gravitational 
dynamics, based on a decomposition relative to a chosen 4- velocity field u^. The 
fundamental tensors are 

where h^i, projects into the instantaneous rest space of comoving observers, and e^^a 
is the projection of the spacetime alternating tensor ri^^ap = — -\/^.9'5°[m^^'^'^^"'^'^^] ■ 
The projected symmetric tracefree (PSTF) parts of vectors and rank-2 tensors are 

^(M> = /i/K , 5(^.> = {ha'^K)^ - ^/i"^ V}^a/3 . (A.2) 

The skew part of a projected rank-2 tensor is spatially dual to the projected vector, 
S^ — j^A"^""^ ' ^^^ then any projected rank-2 tensor has the decomposition 
S^v = lShfj_„ + e^^aS"' + 5'(^^) , where S = Saph"^ . 

The covariant derivative V^ defines 1-1-3 covariant time and spatial derivatives: 

>•..., = 7/"v, j'^-...,, v„ j^-..., = hjh\ ■ ■ ■ K^Vpr-...r.{A.i) 

The projected derivative V^ defines a covariant PSTF divergence, V''V^ , \i"^S^„, and 
a covariant PSTF curl, 

curl V,, = e^.eV'-l/" , curl S,,, = s^p^^V^S,)'^ . (A.4) 

The relative motion of comoving observers is encoded in the PSTF kinematical 
quantities: the volume expansion rate, 4-acceleration, vorticity and shear, given 
respectively by 

e = V^M,„ Af, = u^, uj^ = curlw^, a^„ = ^{t^u^^y (A. 5) 

Thus 

V^w^ = -Qhf_,^ + e^^qw" -I- ai_,;, - A^u„ . (A. 6) 

A key identity (valid in the fully nonlinear case) is 

curl V^^ EEE e,,^„ V'' V" V = -a^-w^ , (A.7) 

which shows that curl grad is nonzero in the presence of vorticity (a purely relativistic 
feature, with no Newtonian analogue). A crucial nonlinear commutation relation for 
scalars is 

V(^-^) ~ ^mV- = ^^M - [l^h^^ + '^p. + e,..K^" j V^V^ . (A.8) 
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The PSTF dynamical quantities describe the sources of the gravitational field: 
the (total) energy density p — T^^.u^u'^, isotropic pressure p — i/i^^T^", momentum 
density 9^ = — T^/^^^u"^, and anisotropic stress tt^^ = T^^ij/), where T^i, is the total 
energy-momentum tensor. The locally free gravitational field, i.e. the part of the 
spacetime curvature not directly determined locally by dynamical sources, is given 
by the Weyl tensor C^^afi- This splits into the PSTF gravito-electric and gravito- 
magnetic fields 

^^. - C^a.0«"w^ , H^, = ]^e^^pC''f',^u\ (A.9) 

which provide a covariant description of tidal forces and gravitational radiation. 
The Ricci and Bianehi identities, 

V[^V,,]Ma = R^,uafiu'\ ^'^C^^^afi = -V[^|i?^]„ - -Rgu^aJ-, (A.IO) 

produce the fundamental evolution and constraint equations governing the covariant 
quantities. Einstein's equations are incorporated via the algebraic replacement of the 
Ricci tensor 

i?^"^ = T^"" - -Ta,°'g'"', (A.ll) 

where T'^'^ is the total energy-momentum tensor, including a cosmological constant 
component, Tj^" = — Ag^'^. 

The resulting equations, in fully nonlinear form and for a general source of the 
gravitational field, are: 
Evolution: 

p+{p + p)e + V^q^ = -2A^q^ - a^-'Tr^, , (A.12) 



e + ^e2 + -(p + 3p)-v^A^ 



i„o 1 

^-<J^,a'"' + 2u;^u;^ + A^A'', (A.13) 

4 -. . 

= -^^.9" + e^.aw V - ^%. , (A.14) 

2 1 

oj{^,) + -Qi^i^ + -cuT\Af^ = a^tyUj'' , (A. 15) 

2 1 ~ 

= -cr«(M^i>>" - ^(p^^u) + A(^,A^) , (A.16) 



1 



E{^,u) + 9^A"^ - curl H^,„ + - (p + p)ct 



H{fj,v) + QH^i, + curlE"^^ - -eurl-TTp^ 

- 2A"e^fi(^E,f - ^uoi^^^q,,) + ^(T"(^£,),^g^ . (A.18) 



Inhomogeneity and the foundations of concordance cosmology 
Constraint: 



V^LU,, = A^^^ , 



2 ~ 



1^„ 1^ 1 

2 



v'E^, + -v^^^, - -v^p + -eq^ 



e^.aCT'-^i/^'^ - m.^^Uj" + id,,,^'^ - ^e^.a^'g" : 



'^"H^y + -curlq^ - (p + p)w^ 



-Sfj.iyaCr'^flE' 



" „ T^o^P 






SEni^Uj'^ 



-.T^Liv^ 
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(A.19) 
(A.20) 
(A.21) 

(A.22) 
(A.23) 



The energy and momentum conservation equations are the evolution equa- 
tions (A.12I and (A.14|. The dynamical quantities /c,p, (?^,7r^i/ in the evolution and 
constraint equations ( A.12|-( A.23) are the total quantities, with contributions from 
all dynamically significant particle species. Thus 



(A.24) 



Tr=p}n^n^+p}hr + 2q* 



iP-,.^) 



T '^ 



(A.25) 



where I — r,n, b, c, A labels the species. The asterisk on the dynamical quantities 
p*j,--- is intended to emphasize that these quantities are measured, not in the u^^- 
frame, but in the /-frame, whose 4-velocity is given by 

< = 7/ ("^ + <) , <"m = , 7/ = (1 - viy'^' ■ (A.26) 

The fully nonlinear equations for the / dynamical quantities as measured in the 
fundamental u^-frame are: 

PI =P*i + {iWi {P*i + P*i) + ^liqT^i^ + 7r*rvi^vi,} , (A.27) 

PI =P*i + T^^llvfip} +p}) + 2jiq*/vi^ + Tr*j'"'vif,vi^j , (A.28) 

9/ = qT + ip*i+p*i)vj + |(7/ - 1)9/'' - iiqT'^iv'^'' + iiqT^ii^^l: 

+ ih] {P*i+P*i) vt + T^*r^i- - ^*r^i.vio.u^'} , (A.29) 

- 2-fiqj°'viaU^'^vy - -TT*j"^viaVifih'^'' + 7j {p*j +p*j) vfvyj . 

(A.30) 

The terms in braces are the nonlinear corrections that vanish in the standard perturbed 
FLRW case. The total dynamical quantities in Eqs. ( A.12[ )- |A.23[ ), are given by 

p-J2pi^p-J2p'^ i'-J2ii^ ^"' = J2<- (a-31) 

III I 

Assuming that the species are non-interacting, they each separately obey the energy 
and momentum conservation equations (A. 12) and (A. 14): 

Pi + {pi +Pi)& + V^g^ = -2A^,gj - cr^i.Tr^'' , 



(A.32) 
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4""^ + ^eq^ + {pi + P/)^^ + v'^p/ + v.^r 

= -fT^gJ' + £%„L.''g? - A^TT^ , (A.33) 



where the /-quantities are given by Eqs. ( A.27)-(A.30). 

The Ehlers-EUis covariant kinetic theory description starts by sphtting the photon 
4-nionientuni as 

pA' ^ E{u^' + e^) , e^e^ = 1 , e^'u^ = . (A.34) 

Here E — —u^p^ is the energy and e^ — p^^^ /E is the direction, as measured by a 
comoving fundamental observer. Then the photon distribution function is decomposed 
into covariant harmonics via the expansion 

f{x,p) = f{x,E,e)^F + F^e'' + F^.e^'e" + ■■■ 

= Y^FMA^.E)e^'''-\ (A.35) 

where Mi = IJ-ifJ-2' ■ ■ l^e and e^^* = e^^ • • • e^^ The multipoles Fm^ are a covariant 
alternative to the usual expansion in spherical harmonics. They are PSTF: 

F^..., = F(^...,) ^ F^..., = F(^...,),F^...,u'' = = F^...^ph"^. (A.36) 

The first 3 multipoles determine the radiation energy-momentum tensor, 

rrix)^ fp^Yf{x,p)d'p 



= pru'^u'' + Iprh^" + 2qi''u''^ + Tl^" , (A.37) 

where d^p = EdEd^ is the covariant volume element on the future light cone at event 
X. It follows that the dynamical quantities of the radiation (in the u^-frame) are: 

Pr 



47r / E^F dE, qi^^ — E^ F^" dE , 

Jo 3 Jo 



Q /"OO 

<"' = —/ E^F^"' dE . (A.38) 

We extend these dynamical quantities to all multipole orders by defining the brightness 
multipoles 



U^,...^,= E'F^,...^,dE, (A.39) 

so that 

1 3 1 "i 

n = -p., n- = -,^n- = -.r. (a.40) 

The coUisionless Boltzmann (or Liouville) equation is 

^.,.^-r...„V^ = o, (A.«) 

where p^ = dx'^/du, and we neglect polarization. The covariant multipoles of df /dv 
are given by 



Inhomogeneity and the foundations of concordance cosmology 



24 






[E'+^'F, 



vMt\ 



nv 



2£ 



(2£ + 3)^ 



p-\l2 p-Al^p 



v{Mi 



''m.) 



(M<._ 



/i«-iM^; ' 



(A.42) 

where a prime denotes d/dE . This is a fully nonlinear expression. 

Multiplying Eq. (A.42) by E^ and integrating over all energies leads to the 
brightness multipole evolution equations: 

4 
3 



= n^Mf) + o ®nM£ + V(^^nMf_i) + 



V i-i-uMi 



(2f + 3) 



^W £i/a(^fnA,fj_i) 



(2^ + 3) 

, (^-l)(^+l)(^ + 2) 
(2£ + 3)(2^ + 5) 



' + ^^^^ ^U'^n.M, + (^ + 3)A<^,nM,_,) 



J-i/aMf 



5f 



-CT''(^,nM,_i)i. - (£ + 2)CT(^^^^_jnM,_2> 



(A.43) 



(2£ + 3) 

Once again, this is a fully nonlinear result. The monopole evolution equation is just 
the energy conservation equation, i.e., Eq. (A.32) with I = r, and the dipole evolution 
equation is the momentum conservation equation (A. 33 1, with I — r. The quadrupole 
evolution is given by Eq. Q. 
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